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Abstract 

Based on density functional calculations, optimized structures of graphite oxide are found for various 
coverage by oxygen and hydroxyl groups, as well as their ratio corresponding to the minimum of total energy. 
The model proposed describes well known experimental results. In particular, it explains why it is so difficult 
to reduce the graphite oxide up to pure graphene. Evolution of the electronic structure of graphite oxide 

OO . with the coverage change is investigated. 
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^ ■ 1 Introduction 

^ , Despite the fact that graphite oxide (GO) was first derived more than a century ago^'^ its structure and 

chemical composition remains not quite clear yet. GO can be used for production of graphite nanoparticles 
and an insulating material for nanodevices^'^. Recently, after discovery of extraordinary electronic properties of 
single-layer carbon, graphene^"®, and successful exfoliation of layers in GO^'"'^" it is considered as a perspective 
source of a "cheap graphene""^. Direct structural information about GO can be hardly obtained (for the most 
of structural methods the use of bulk crystals are desirable whereas the GO exists mainly in solutions) which 
C/3 , makes especially important theoretical modeling of its structure and properties. 

Original methods of preparation of GO have been modified afterwards^ ^"^"^ which allows to vary a bit 

its chemical composition. In Nakajima and Matsuo paper^^ the chemical compositions of GO derived by 

the methods developed by Brodic"'^ and Staudenmaier^ were determined as C8H2.54O3.91 and C8H4.61O6.70, 

respectively. This means that both hydroxyl groups connected with single carbon atoms and oxygen atoms 

C^ ' connected with two carbon atoms present in GO (see Fig. [ij. According to Szabo and co-workers^^, the chemical 

composition of different samples of GO varies from C8H1.20O3.12 to C8H1.60O3.92, according to Hontora-Lucas 

et al. paper^^ - from C8(OH)i.380o.63 to C8(OH)i.640o.79, and, according to Cassagneau and co-workers^^ - 

from C12HO2 to C15H3O4. In general, one can conclude that the chemical composition of GO, not considering 

Q ' groups coupled with graphene edges, varies in a range from C8H2O3 to C8H4O5. It is a common opinion (see 

O , the papers cited above) that oxygen in GO mainly presents in hydroxyl groups or in bridges connecting two 

carbon atoms in graphene layers whereas the amount of carboxyl, as well as carbonyl, groups is relatively small. 

^vq , Thus, the two limiting compositions of GO can be presented as C8(OH)2 and C8(OII)40. All chemical formulas 

^ ' of GO obtained experimentally manifest the coverage of graphene between 25% and 75%, which means that, 

■"nJ" , at least, quarter of carbon-carbon bonds in graphene layer are double bonds whereas the rest is single bonds 

2P ■ like in diamond. Indeed, both XPS^'^'^^ and infrared^^ spectroscopy data confirm coexistence of sp^ and sp^ 

electron configurations of carbon. 

Based on experimental data^^ and additional measurement results, a model of GO has been suggested by 
^^ I Nakajima and co-workers"'^*. According to this model, the GO structure is intermediate between two ideal 

structures, C8O2 and C8(OH)4 (Fig. [^t and h respectively). Later models^^'^° differ mainly by assumptions 
r^ ] concerning edge groups. As a result, the GO structure is described as a combination of completely covered and 

completely uncovered stripes of graphene which is confirmed also by recent theoretical results^ "'^'^^. Recently, 
mechanical properties of GO have been simulated, based on modeling of nanoobjects functionalized by oxygen 
^^ ■ from one side^^ or on experimental data on chemical composition of GO^"'. 

However, due to the stripes of uncovered graphene, GO should be conducting, according to these models. At 
the same time, experimentally GO becomes conducting only after a very strong reduction'*'^^ whereas typically 
GO is insulating. Here, based on density functional calculations, we formulate a model of insulating GO. We 
investigate also a transition to conducting state at the reduction and explain why it is so difficult to clean GO 
completely and to derive pure graphene from GO. 

2 Computational Method 

Some general factors determining chemical functionalization of graphene have been investigated in our previous 
work^^ using hydrogenization as an example. First, graphene is a very flexible material, and a chemisorption 
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of even single hydrogen atom leads to essential distortions of the graphene sheet with a radius approximately 
5 A, and these lattice distortions are of crucial importance for energetics of the process. Second, for the 
chemisorption of two hydrogen atoms the configuration where they are bonded with two neighboring carbon 
atoms from opposite sides of the sheet turns out to be the most energetically favorable. Third, complete coverage 
by hydrogen provides the global minimum of energy. These features are relevant, as we will show here, also for 
the case of GO. 

We used the pseudopotential density functional SIESTA package for electronic structure calculations^^'^^ 
with the generalized gradient approximation for the density functional^^ , with energy mesh cutoff 400 Ry, and 
fc-point 11x11x1 mesh in Monkhorst-Park scheme^". During the optimization, the electronic ground states 
was found self-consistently by using norm-conserving pseudopotentials to cores and a double-C plus polarization 
basis of localized orbitals for carbon and oxygen, and double-C one for hydrogen. Optimization of the bond 
lengths and total energies was performed with an accuracy 0.04 eV /A and 1 meV, respectively. This method 
is frequently used for computations of electronic structure of graphene^^''^^"'^^. 

To compute the properties of layered GO we have carried out calculations of the corresponding structures 
for the case of 25% coverage by hydroxyl groups. Instead of GGA, we use here the LDA approximation which 
is known to be more accurate to describe interlayer coupling in graphite and other van der Waals systems^^. 
The basis for carbon atoms was optimized as proposed earlier for pure graphite^^. 

Chemisorption energies were calculated by standard formulas used, e.g., earlier for the case of chemisorption 
of hydrogen on graphene^^ and solution of carbon in 7-iron'^^. Thus, energy of chemisorption of single oxygen 
atom at eight carbon atoms (Fig. 2a) is calculated as E/orm=Ec80-Ec8-Eo2/2 where EcgO is the total energy 
of the supercell found by self-consistent calculations after optimization of geometric structure, E^g is the total 
energy of graphene supercell, and E02 is the energy of oxygen molecule. For the case of hydroxyl group, instead 
of oxygen, its energy was calculated with respect to water in gaseous phase: Yioh=^H20-^H2I'^- Alternatively, 
the chemisorption energy can be calculated as Eoff=E//2o/2+Eo2/4. These two expressions estimate the 
chemisorption energy from above and from below (see Fig. [3}d where the results corresponding to the first 
and to the second expression are shown by dashed green and dotted blue, respectively). To be specific, in 
further discussions we will use the first estimation (the dashed green line). Actually, the chemisorption energy 
for GO containing both oxygen and hydroxyl groups depend on its chemical composition. For example, the 
chemisorption energy of oxygen and OH group in the system C8(OH)40 are ^chem = ^Cfi{OH)4, + E02/2 - 
Ec8(OH)40 and ^chcn = {^Cs(oh)20 + 2Eoh - Ece(off)4o)/2, respectively. 

To check an accuracy of the method used we have calculated formation energy of the water from molecular 
oxygen and hydrogen in gaseous phase. We have found the value 213.4 kJ/mol which is rather close to the 
experimental value 241.8 kJ/mol. Underestimation of the energy by approximately 10% is typical for GGA 
calculations'^^. Also, we have calculated equilibrium interatomic distances for graphene, molecular oxygen, 
hydrogen, and water, as well as interlayer distances in graphene. When drawing the pictures of density of 
states, a smearing by 0.3 eV was used. 

3 Results and Discussion 

We start our simulations with the case of oxygen chemisorption, then consider the chemisorption of hydroxyl 
groups and, at last, investigate their various combinations. Let us consider first a supercell of graphene contain- 
ing 8 carbon atoms, the chemisorption of two of them corresponding 25% coverage. In contrast with hydrogen, 
oxygen forms a bridge between two carbon atoms, as shown in Fig. [T^. As well as for the case of hydrogen, 
the chemisorption leads to distortions of graphene sheets when the atoms coupled to oxygen are shifted up 
and their neighbors are moved in the opposite direction. This makes chemisorption of the next oxygen atom 
from the opposite side of graphene sheet the most energetically profitable (Fig. [T^). Various oxygen configura- 
tions for various coverage are sketched in Fig. [2K-e and the computational results for carbon-carbon distances, 
chemisorption energies, and electron energy gaps are presented in Fig. [31 One can see from Fig. [3^ that the 
length of bonds between functionalized carbon atoms grows from the standard value for graphene, 1.42 A, to 
the standard value for diamond, 1.54 A, at the coverage increase which corresponds to the transition from sp^ 
to sp^ hybridization of carbon atoms. The chemisorption energy increases in absolute value with the coverage 
increase and the most stable is the configuration displayed in Fig. [T^. The gap in electron energy spectrum 
opens starting form 75% coverage where its value is 1.8 eV; with the coverage increase, it grows up to 2.9 eV. 
Hydroxyl groups are bonded with graphene similar to hydrogen^^, that is, they sit at neighboring carbon 
atoms from opposite sides of the graphene sheet (Fig. [TJd). Distortions of the sheet is stronger in this case 
than in the case of hydrogen, partially, due to interaction between the hydroxyl groups leading to ordering 



of distortion (see Fig. [Ib,c). Various calculated configurations are sketched in Fig. Ef-g- The chemisorption 
energy was calculated with respect to water which is probably the most informative to consider reduction of 
GO. With the coverage increase, the carbon-carbon distance growth, as shown in Fig. [3^. For the case of 
complete coverage it turns out to be larger than a standard one for a single bond {sp^ hybridization) which 
means a situation close to break of the graphene sheet. In contrast with the cases of hydrogen and oxygen, the 
chemisorption energy is not monotonous as a function of coverage reaching the minimum at 75% which should 
correspond, therefore, to the most optimal configuration (Fig. [Jd). 

Let us consider now a general case of functionalization of graphene by both hydroxyl groups and oxygen 
atoms. Typical combinations are shown in Fig. U] Total energy calculations demonstrate that for all the 
combinations under consideration the chemisorption energy per hydroxyl group is 60 meV lower, and per 
oxygen atom 30 meV lower than for the pure cases with the same degrees of coverage. Thus, mixed coverage is 
energetically favorable in these cases (coverage between 25% and 75%) diminishing the energy of both O and 
OH groups. For coverage less than 25% the chemisorption energy for hydroxyl groups turns out to be lower 
than for oxygen, therefore one can conclude that GO with 25% coverage contains only OH groups whereas some 
oxygen atoms can appear only as edge groups. Optimal configurations for 25%, 50%, and 100% coverage are 
shown in Figs. H^, e and c, and h, respectively. One should stress that a structure with staggered stripes of sp^ 
and sp^ carbon atoms is formed at 50%, in agreement with the previous works^^'^^. For the maximal coverage, 
as well as for the case of OH groups only, carbon-carbon distances exceed 1.54 A which makes, again, 100% 
coverage less favorable than 75% one. The most stable configuration of GO is presented in Fig. [TJ:. 

As a result, one can suggest the following chemical formulas for GO with various coverage: 25% - C8(OH)2, 
50% - C8(OH)20, and 75% - C8(OH)40. They are rather close to the formulas suggested by experimentalists 
and discussed in the Introduction. Minor discrepancies can be related with the presence of some small amount 
of carboxyl and carbonyl groups, as well as atomic oxygen adsorbed at the edges of GO, as was discussed in 
detail within the model proposed by Lerf et al^^'^*'. 

Electron densities of states for GO are presented in Fig. [5) One can sec that the energy gap varies between 2.8 
eV and 1.8 eV at the decrease of coverage from 75% to 50%. At further reduction of GO it becomes conducting, 
according to our calculations. It seems to be in agreement with the available experimental data'''^^'^^. 

The chemisorption energy difference per group for 25% and 75% coverage is less than 1 eV (Fig. [3^) which 
explains a possibility of partial reduction of GO, both thermally and chemically. Actually, the carbon to oxygen 
ratio 4:1 considered above is a bit larger than experimental values for strongly reduced GO^^'^^''^*''^^ and almost 
twice larger than the maximal ratio 10:1^^. To study dependence of the chemisorption energy on the C:0 ratio 
we have performed calculations for the cases of two hydroxyl groups (see Fig. [2f) per supercells containing 
8, 18, 32, 50, and 72 carbon atoms, the latter case corresponding to the C:0 ratio 32:1. The computational 
results are presented in Fig. [6^. One can see that the chemisorption energy is weakly concentration-dependent 
between the ratio values 4:1 and 16:1 whereas for smaller concentrations of hydroxyl groups it decreases roughly 
twice, between 16:1 and 25:1. A weakening of chemical bonding between OH groups and graphene at small 
concentrations (C:0 ratio from 25:1 in comparison with 16:1) is connected with essential changes of the electronic 
structure. At very small concentrations, the latter is more similar to that of pure graphene (see Fig. [Dd). It can 
be caused by the changes of distances between OH groups which is 17 A for the carbon-to-oxygen ratio 16:1. 
It was shown at the simulation of hydrogenization of graphene^^ that typical radius of interaction between the 
hydrogen atoms is about 8 A, and the defects can be considered as independent ones for larger distances. For 
the case of OH groups, the distortions of graphene sheet is larger than for the case of hydrogen and therefore 
interaction between them is still essential for the ration 16:1 whereas for smaller concentrations the hydroxyl 
groups can be considered as almost non-interacting. In real experimental situation where the ratio 10:1 has 
been reached^^ finite-size effects of the GO flakes can be important. Indeed, the size of these flakes is smaller 
than for the case of graphene^'^^ and various groups can be chemisorbed at the edges. Also, the flakes can 
contain various topological defects^°'^'^ which can also change local chemisorption energy. 

At last, let us discuss the cases of bilayer and periodic (graphite-like) GO. To this aim we have carried out 
calculations for corresponding structures with 25% coverage by hydroxyl groups. We consider the Bernal (AB) 
stacking, similar to pure graphite, which was observed also in GO"'^^. The optimized structure is shown in Fig. 
[71 The width of the layer was found to be, in both cases, about 7 A, as well as interlayer distances, which seems 
to be in a good agreement with experimental data^^'^^''*^. To calculate interlayer coupling energy per carbon 
atom we have computed the energy differences between single layer and periodic structure. For the case of pure 
graphite, it equals 32 meV, in a good agreement with the experimental value 35 meV*^. For the periodic GO 
structure and for GO bilayer, the corresponding values turn out to be 17 meV and 6 meV, respectively. This 
decrease of the energy explains possibility of easy exfoliation in GO^'^*^. Due to weakness of interlayer coupling, 
the electronic structure of GO is almost identical for single layer, bilayer, and periodic structure, in contrast 



with the cases of pure graphene^^. 

4 Conclusions 

To summarize, we have proposed a model structure of GO which seems to be consistent with all known experi- 
mental data. We have demonstrated, in particular, that (i) 100% coverage of GO is less energetically favorable 
than 75% (ii) functionalization by both oxygen and hydroxyl groups is more favorable for coverage than 25% 
than by hydroxyl groups only (iii) a reduction of GO from 75% to 6.25% (C:0 ratio 16:1) coverage is relatively 
easy but further reduction seems to be rather difficult, and (iv) GO becomes conducting at the coverage 25%, 
being an insulator for larger coverage. 
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Figure 1: The most stable configurations of grapliene functionalized by oxygen only (a), hydroxyl groups only (b), 
and both oxygen and hydroxyl groups (c). Carbon, oxygen and hydrogen atoms are shown in green, blue, and violet, 
respectively. 
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Figure 2: A sketch of functionalization of graphene by (a-e) oxygens (blue circles), and (f-g) hydroxy groups (green 
circles) . 
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Figure 3: Dependence of carbon-carbon bond length (a), chemisorption energy (b), and electron energy gap (c) on 
coverage (see the text). 
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Figure 4: A sketch of functionalization of graphene by oxygens (blue circles) and hydroxyl groups (green circles). 
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Figure 5: Electronic densities of states for the most stable configurations at various degrees of coverage. Number of 
hydroxyl groups and oxygen atoms per Cg is shown. 
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Figure 6: (a) Chemisorption energy of OH groups as a function of C:0 ratios; (b) Total densities of states per atom 
for C:0 ratios 16:1 (solid red line) and 32:1 (dashed green line). 
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Figure 7: Optimized geometric structures of strongly reduced GO. Numbers are distances, in A, for the periodic 
structure (and for bilayer in parentheses). Right upper corner: a top view. Carbon, oxygen and hydrogen atoms are 
shown in green, blue, and violet, respectively. 
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